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Abstract 

EGF and its receptor EGFR serve as a paradigm for signaling In cell, molecular and tumor biology. EGFR Inhibitors, drugs 
targeting the intracellular kinase activity and antibodies targeting the extracellular llgand binding, are used to treat breast, 
lung, colon and other cancers. Nominally affecting the same target, Inhibitors have different effects, suggesting that use of 
inhibitor combinations may provide beneficial In cancer treatment. To explore the specific and the common transcriptional 
effects of EGFR inhibitors, we present metaanalysis of 20 individual studies comprising 346 microarrays. We identified 
specific gene subsets regulated by kinase Inhibitors, those regulated using antibodies and by suppressing EGFR expression 
using miR-7. Unreported before, the Inhibitors prominently induce lysosome components. All Inhibitors rely on related sets 
of transcription factors and protein kinases, both for transcriptional Induction and suppression. However, we find that 
Gefltlnib suppresses apoptosis inhibitors, while Inducing cell-cycle inhibitors; conversely, Eriotinib suppresses cell-cycle and 
cell migration genes, while inducing proapoptotic genes. EGFR-targetIng antibodies specifically suppress cell motility, 
developmental and differentiation processes, while inducing the contractile apparatus. mlR-7, distinctively, suppresses cell- 
cycle genes, while inducing transcription machinery. These metaanalysis results suggest that different inhibitors have 
overlapping but quite distinct effects In target cells. Judicial use of EGFR-targetIng combinations, i.e., simultaneous use of 
antibodies and multiple kinase Inhibitors, may provide more effective cancer treatments with fewer side-effects and avoid 
development of resistance. We expect, moreover, that specific drug combination treatments can be fine-tuned to achieve 
specific, personalized results. 
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Introduction 

Epidermal growth factor, EGF, afTects almost all cell types, 
including eponymous epidermis; its signaling is deregulated in 
many pathological conditions [1]. EGF and its receptor EGFR 
constitute, arguably, the most studied model of cellular signaling 
[2-4]. EGFR responsive signaling pathways include GRB2/ 
MAPK, PI3K/AKT, STATs, PLC/PKC, and transcription 
factors API, Myc, Egrl etc. [1,2,4—6]. The EGF-regulated genes 
promote ceU-cycle and proliferation, protein synthesis, migration, 
adhesion, ECM remodeling, angiogenesis, tumorigenesis and 
metastasis; conversely, apoptosis and terminal difTerentiation are 
usually inhibited [7-9] . 

Activated EGFR typifies numerous epithelial malignancies, 
including cancers of the breast, lung, colon, head-and-neck, 
pancreas etc. [7]. Therapies that inhibit EGFR became a 
paradigm for targeted treatment of human cancers and use 
inhibitors of EGFR kinase, Gefitinib and Eriotinib (a.k.a. Iressa 
and Tarceva, resp.), or antibodies Lapatinib, Cetuximab, 
Panitumumab, Zalutumumab, Niniotuzumab and Matuzumab 
[7,10]. They can induce tumor regression avoiding some adverse 
effects of chemotherapy. Drawbacks of EGFR inhibitor therapies 
are cardiac and renal side-efTects, skin toxicity, and intrinsic or 
acquired resistance to therapy; these limit the duration or dosage 
of treatment [,5, 11]. 



Whereas all these agents target the same protein, EGFR, 
diflferent inhibitors use different mechanisms and have different 
effects [12]. For example, Gefitinib and Eriotinib compete with 
ATP and inhibit receptor autophosphorylation, retaining effec- 
tiveness against constitutively active kinase mutants. Antibodies 
bind the extracellular domain of receptor, occluding ligand 
binding, preventing receptor dimerization and activating host 
immune responses [12,13]. 

Many studies used transcriptional profiling to define cellular 
responses of targeting EGFR. However, the use of different agents, 
microarray platforms and experimental protocols makes it difihcult 
to characterize the commonalities and the particulars of EGFR 
inhibition. Our objective here is to use metaanalysis for a 
comprehensive investigation of transcriptional data. We metaana- 
lysed 20 published transcriptional studies, comprising 346 
microarrays, using free, readily available computer programs, 
e.g., RankProd [14]. We determined the ontological categories 
overrepresented in the regulated genes and identified potential 
protein kinases and transcription factors involved. 

The results describe large lists of over 2537 suppressed genes 
and 2251 induced by EGFR inhibitors, with high statistical 
significance. They identify crucial differences in the genes 
regulated by antibodies and by kinase inhibitors and specifically 
the consequences of Gefitinib vs. Eriotinib treatments. We also 
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demonstrate the great advantage of metaanalysis over single 
studies. The work can serve as a paradigm for integration and 
metaanalysis of transcriptional data in public repositories. 

Methods 

To identify transcriptional studies of EGFR inhibition, we 
searched, PubMed GEO_DataSets [http://\™'w.ncbi.nlm.nih. 
gov/pubmed] and EMBI-EBI_ArrayExpress [http://www.ebi.ac. 
uk/arrayexpress], and flagged 50 and 44 data sets, respectively. 
The results were individually screened yielding 20 studies 
(Table 1). In 2 of these, GSE17948 and GSE40130, different 
inhibitors types were used, which were analyzed separately. Thus, 
we analyzed 22 different sets that compared direcdy EGFR 
inhibitor-treated vs. untreated samples. These are listed in Table 1 
and described in more detail in Supplement Text SI. Two 
additional studies were identified: GSE6128 compared inhibitor- 
treated cells with the same cells after rcm()\ al of the inhibitor, and 
GSE 17498 used a platform with fewer gene targets; these two 
studies were not used in our metaanalysis. Most datasets used 
Affymetrix microarray platforms (1 1 studies), lUumina and Agilent 
microarrays were used in 6 and 3 studies, respectively (Table 1). In 
several studies multiple cell lines or tissues were treated, in these 
cases each cell line was analyzed separately. 

We recently published detailed metaanalysis protocols for 
finding relevant studies in public repositories, downloading data 
files, microarray quality control, cross-referencing and merging 
different platforms, dealing with empty cells and negative 
numbers, selecting data subsets and differentially expressed genes 
using RankProd, gene lists annotation using the Lists2Networks 
algorithm, clustering and comparing lists and promoter analysis 
[15]. The workflow diagram describing meta-profiling and cluster 
analysis is shown in Fig. 1 Briefly, for the data deposited in studies 
using the AffymetrLx arrays we downloaded CEL files and then 
processed these using RAIAExpress software [16]. To combine 
data from different platforms we used a non-parametric approach, 
RankProd, which identifies differentially expressed genes, both up- 
or down- regulated, based on the estimated percentage of false 
predictions [14]. Lists of official symbols of regulated genes was 
submitted to the Lists2Networks analysis program [17]. The 
program compares lists for mutual overlaps within ontological 
categories, kinase targets, and disease-associated genes, providing 
statistical evaluation of the overlaps. For annotation we used the 
DAVID program suite [18]. 

Results 

Identification and characterization of data sets 

Because EGFR activation often plays important role in cancers, 
many teams used various inhibitor agents, microarray platforms 
and experimental approaches in transcriptional studies of EGFR 
inhibition. We identified such studies, and then, used RankProd to 
combine data for metaanalysis. Several studies used multiple cell 
types, which were compared separately, producing 43 individual 
pair-wise comparisons (Table 1). 

One of strengths of the metaanalysis approach is its flexibility, 
which allows direct comparisons of subsets of studies. For example, 
we selected and separately identified regulated genes in studies 
using inhibitors of EGFR kinase, ones using antibodies, and one 
study using miR-7 to suppress EGFR production. The numbers of 
differentially expressed genes in study subsets are summarized in 
Fig. 2. 

Using DAVID [18], we explored the ontological categories and 
signaling pathways. The complete annotated hst of genes induced 



and suppressed with p-value better than 10 * from all 346 samples 
is given in Tables SI, S2. In addition to providing annotations, 
DAVID estimates which ontological categories and pathways are 
over-represented in the submitted lists. In Table S3 we provide the 
top 100 categories in each list, including all 346 samples, and 
subsets for specific inhibitors. Because DAVID scans many 
different data bases to identify over-represented categories, the 
resulting charts contain many redundancies. For example, "cell- 
cycle", "cell division", "mitosis" and "M-phase" are listed as 
different categories. Such redundancies are largely condensed by 
'clustering', which combines categories tfiat comprise largely 
overlapping sets of genes [18]. Therefore, below we focus our 
analyses on clusters of differentially expressed genes. 

EGFR kinase inhibitors 

Several studies used specific kinase inhibitors that compete for 
binding the ATP pocket of the kinase, Gefitinib, Erlotinib, 
Lapatinib and Tyrpliostin AG1478. While these agents are 
expected to act by similar mechanism, i.e., inhibiting the kinase 
enzymatic activity, subde differences in their effects were seen 
[12]. The doses of the inhibitors varied among studies, depending 
on the inhibitor, target cell type etc., but were generally aimed at 
fuU, saturating doses that completely inhibit the EGFR kinase. We 
identified a subset of studies using all kinase inhibitors (i.e., 
excluding antibodies and miR-7), a subset using Gefitinib, the 
largest such subdivision, one using all inhibitors except Gefitinib, 
and one using Erlotinib. We then identified the commonalities and 
particulars of these subsets and compared them to the complete 
346 samples set. 

The ontological category 'nuclear lumen' is statistically the most 
prominent one in the suppressed genes (Table 2). Protein 
modification and metabolism processes, marked 'P' in Table 2, 
are prominent in all suppressed categories, as are transcription and 
translation, 'T'. Regulators of apoptosis, both positive and 
negative, are suppressed by EGFR inhibitors. Interestingly, 
migration and cell-cycle categories, 'M' and 'C', did not reach 
statistical significance in the Gefitinib-suppressed genes, although 
they are very prominent in the Erlotinib-suppressed genes 
(Table 2). 

The most prominent categorv' induced by EGFR kinase 
inhibitors is lysosome, 'L', the (^toplasmic membrane-bound 
organelle containing hydrolytic enzymes of intracellular degrada- 
tion (Table 3). Inhibitors of transcription and translation are 
induced, as are positive regulators of apoptosis; these effects neatiy 
dovetail with their opposites in the suppressed categories. 
Interestingly, the cell-cycle inhil)itors are induced by Gefitinib: 
apparently, the mechanisms of shutting down the cell-cycle differ 
among EGFR kinase inhibitors: while Erlotinib suppresses cell- 
cycle proteins, Gefitinib induces cell-cycle inhibitors. This 
interesting observation deserves additional studies because, if 
confirmed, it could have important implications for cancer 
treatments. 

Steroid hormone receptor activity was conspicuous, particularly 
in the Gefitinib-induced set (Table 3). While its significance is 
currently unknown, we note the nexus between EGFR inhibition 
and androgen activity [19], that anti-estrogens are used combined 
with Gefitinib against breast and lung cancers [20], whereas 
corticosteroids are used to treat side-effects of Gefitinib [21]. 

Ontological category 'neural tube formation', 'E' in Table 3, 
was also prominent; PubMed search of 'Neural tube AND 
Gefitinib' gave no results. 
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Table 1. Data sets used in the metaanalysis. 



Data set 


Platform 


Set 


Control+treated 


Target cells, 
pretreatment 


Inhibitor 


1 GSE11729 


Affy HG U133 Plus 2 


1 


4-H2 


HI 299 


gefitinib 






2 


12-1-12 


HI 299, EGF 


gefitinib 






3 


2+6 


L858R 


gefitinib 






4 


6+6 


L858R, EGF 


gefitinib 


2 GSE16179 


Affy HG U133 Plus 2 


5 


3+3 


BT474 


lapatinib 






6 


3+3 


BT474-J4, 


lapatinib 






7 


3+3 


BT474-J4, forteinib 


lapatinib 


3 GSE20854 


Affy HG U133 Plus 2 


8 


4+4 


HecSOco 


gefiitinib 






9 


4+4 


Ishikawa H 


gefiitinib 


4 GSE2320e 


Affy HG_U133_P!us_2 


10 


1+1 


H322C 


gefiitinib 


5 GSEBOSne 


Affy HG_U133_Plus_2 


11 


4+9 


BT20 


eriotinib 






12 


3+3 


MD 


eriotinib 






13 


3+3 


MCF 


eriotinib 


6 GSE32975 


Affy HG_U133_Plus_2 


14 


2+2 


HaCat, EGF 


gefiitinib 






15 


2+2 


HaCat, serum 


gefiitinib 


7 GSE33658 


Affy HG_U133_Plus_2 


16 


5+6 


breast cancer biopsies 
+/- Gefitinib 


gefiitinib 


8 GSE6521 


Affy HG_U133_Plus_2 


17 


2+7 


MCF7 


AG 1478 






18 


7+7 


MCF7, heregulin 


AG 1478 


9 GSE8141 


Affy HG_U133_Plus_2 


19 


4+4 


MCF7 


gefiitinib 






20 


4+4 


MCF7/HER2-18 


gefiitinib 






21 


4+4 


MCF7/HER2-18, Tamoxifen 


gefiitinib 


10 GSE19S00 


Affy HG_U133A 


22 


4+4 


bronchial epithelial cells 


AG 1478 






23 


3+3 


bronchial e.c, high 
barometric pressure 


AG 1478 


11 GSE17948 


Affy HGU95V2 


24 


3+4 


HMCV-infected monocytes 


AG 1478 






25 


3-Ht 


HMCV-infected monocytes 


antiEGFR 


12 GSE23428 


Agilent custom array 


26 


5+5 


triple-negative breast 
cancer biopsies 


cetuximab 


13 GSE32333 


Agilent-027114 


27 


3-14 


xenografts of A431 


nimotuzumab 






28 


4+3 


xenografts of A431, 
Rapamicin 


nimotuzumab 


14 GSE38302 


Agilent-028004 


29 


l-H 


PC-9 


gefiitinib 






30 


3+3 


PC-9GR 


gefiitinib 


15 GSE38310 


lllumina HumanHT-12 v3.0 


31 


3+3 


HCC827 


eriotinib 






32 


3+3 


HCC827-ER3 


eriotinib 






33 


3+3 


HCC827-T1 5-2 


eriotinib 


16 GSE38376 


lllumina HumanHT-12 v3.0 


34 


3+6 


SKBR3 


lapatinib 






35 


3+6 


SKBR3-R 


lapatinib 


17 GSE40130 


lllumina HumanHT-12 V4.0 


36 


3+3 


FaDu, eriotinib 


miR-7 




lllumina HumanWG-6 v3.0 


37 


2+2 


HNS, eriotinib 


miR-7 


18 GSE19043 


lllumina HumanRef-8 v2 


38 


3+3 


DiFi 


gefitinib 






39 


3+3 


GTL-16 


gefitinib 


19 GSE23175 


lllumina HumanRef-8 v2 


40 


1+15 


MCF10-A 


gefiitinib 






41 


1+15 


MCF10HER 


gefiitinib 


20 GSE34557 


lllumina HumanWG-6 v3.0 


42 


3+3 


human epidermal 
keratinocytes 


AG 1478 






43 


3+3 


human epidermal 
keratinocytes, BMP2 


AG 1478 


Only samples directly comparing EGFR inhibitor-treated with their control samples were used. The comparisons column designates specific sets of treated and control 
samples. Note that GSE40130 used two different inhibitors, which were analyzed separately, 'a' marks the studies used in the 'Not-Gefitinib' set. 
doi:l 0.1 371 /journal.pone.Ol 02466.t001 
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Search GEO datasets and 
ArrayExpress datasets 



Manually select studies 
that directly compare 
EGFR-inhibitedand not 
inhibited samples 



Download 

Affymetrix CEL files 




RMAExpress 





Download non- 
Affymetrix TXT 

files 



RankProd non-parametric 
selection of differentially 
expressed genes 



Annotation and 
ontologies using 
DAVID 



Relations among 

lists using 
Lists2Networks 



Figure 1. Workflow diagram. The steps used in data assembly and 
metaanalysis. 

doi:l 0.1 371/journal.pone.01 02466.g001 



Signaling pathways responsive to l<inase inhibitors 

The kinases signaling downstream from EGFR, especially the 
MAPK3K family and of course, the EGFR itself, are very 
prominent, in both suppressed and induced sets (Table S4). 
Comparing the Gefitinib- and Erlotinib-associated kinases, we see 
extensive parallels, with the notable exception of GSK3A/ 
GSK3B, present in the Gefitinib, but not Erlotinib list (Table 
S4). GSK3 is a multifunctional protein controlling glycogen 
synthase, JUN and other transcription factors, with a role in the 
WNT and phosphoinositide-3-kinase signaling pathways [22]. 
PubMed search of Gefitinib & GSK3A yielded no hits, although a 
nexus between EGFR and the Wnt pathway has been reported 
[23]. 

We used \'cnn diagrams to compare overlaps among genes lists 
(Fig. 3). Genes regulated by individual kinase inhibitors signifi- 
candy overlap with the genes tagged as regulated by all inhibitors. 
306 genes are induced and 247 are suppressed by every kinase 
inhibitor (Fig. 3). We separately analyzed these genes: all kinase 
inhibitors suppressed nuclear materiel, translational machinery, 
protein kinases, cell migration, ECM binding, and regulators of 
apoptosis; conversely, all inhibitors induced DNA binding 
proteins, particularly inhibitors of transcription (Table 4). 

We analyzed separately the sets of genes responding to 
Gefitinib, Erlotinib, and the non-Gefitinib inhibitors (Fig. 3, 
Table S5). A caveat in this analysis is that differentially reaching 
statistical significance does not necessarily mean differential 
transcriptional responses. Negative regulators of transcription are 
suppressed by Gefitinib (Table S5). This is surprising because in 
parallel to other kinase inhibitors, Gefitinib suppresses positive 
regulators of transcription (Tables 2, 3). We hypothesize that 
Gefitinib induces very specific transcriptional responses, over a 
background of general inhibition of transcription. The transcrip- 
tion/translation machinery is, as expected, suppressed by the 



a) 



Differentially expressed genes 
p-value < 0.0001 



4500 



4000 
3500 ^ 
3000 
2500 

2000 
1500 
1000 
500 
0 




All Kinase Gefitinib Not Get. Erlotinib Antibodies mlR-7 
InhSiitors 

□Suppressed ■ Induced 



Figure 2. iVIetaanalysis identifies genes regulated by EGFR inhibitors. Bar graph showing numbers of genes identified as differentially 
expressed with p-values 10 or better on different platforms, or using different EGFR inhibitors. Note that for miR-7 p-value 10"^ was used as the 
cut-off because only 10 samples were available. 
doi:1 0.1 371/journal.pone.01 02466.g002 
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non-Gefitinib inhibitors, while inhibitors of protein transport, 
lysosomal degradation and transcription are induced (Table S5). 

Unanticipated, Gefitinib induces the cell-cycle machinery 
(Table S5). This is an unexpected response to EGFR inhibition 
and we note that Gefitinib, unlike Erlotinib and other kinase 
inhibitors, does not generally suppress cell-cycle genes (Tables 2, 
3). Confirming the above, the non-Gefitinib kinase inhibitors 
(which include Erlotinib) specifically suppressed the cell-cycle 
machinery (Table S5). This observation reinforces the hypothesis 
that Gefitinib, specifically among EGFR kinase inhibitors, may not 
directly inhibit the cell-cycle. Obviously, such contentions need 
direct lab-bench proof 

The transcriptional changes in response to EGFR inhibition 
reflect, presumably, the changes in the activity of transcription 
factors. We identified the transcription factors with binding sites 
statistically overrepresented in the regxilated genes. In general, 
very similar sets of transcription factors appear activated by 
difiFerent receptors. There is an overlap between the transcription 
factors responsible for the induced genes and for the suppressed 
ones (Table S6). 

Antibody inhibitors targeting EGFR 

Antibodies inhibit EGFR by binding to the extracellular domain 
of EGFR, occluding the ligand binding site and preventing 
receptor dimerization. Treatment with antibodies is plagued by 
resistance from mutations that constitutively activate the kinase. 
To analyze the transcriptional effects caused by EGFR-targeting 
antibodies, we cuUed from the total only GSE32975, GSE 17948, 
GSE23468 and GSE323333, 39 samples in total (Table 1). 
Because the number of studies and microarrays using antibodies 
is smaller, the 'enrichment scores' for the ontological categories are 
lower too. 

Specifically, antibodies suppress genes associated with cell 
migration (Table 5). EGFR signaling is involved in actin 
cytoskeleton reorganization; however, the antibody inhibitors, 
but not kinase inhibitors, specifically induce the cytoskeletal 
contractile elements dominated by myosin genes (marked with 
tildes in Table 5). At present we have no explanation why the 
antibody inhibitors preferentially regulate the contractile appara- 
tus. Antibodies do not suppress cell-cycle or apoptosis associated 
genes, which are prominently suppressed by kinase inhibitors (c.f 
Table 2). 

Importantly, none of the protein kinases or transcription factors 
reached statistical significance in association with antibody- 
regulated genes (Table S4). Comparing overlaps among differen- 
tially regulated genes, we find 255 genes specifically suppressed by 
the antibodies, while 112 (40+72) are suppressed in common with 
the kinase inhibitors (Fig. 3). Focusing on these two sets, we find 
that the ontological categories suppressed only by the antibodies 
are associated with developmental processes, specifically with a) 
epidermis/hair, b) reproduction/ovulation and c) vasculature/ 
blood development (asterisks in Table 5). These results are quite 
unexpected and point, on one hand, to serious side-effects of 
EGFR-targeting therapies and on the other underline the 
importance of EGFR signaling in the homeostasis of these self- 
renewing tissues. 

The suppressed processes in common with the kinase inhibitors 
comprise migration and responses to steroid hormone; the induced 
ones include apoptosis regulation, but, importantly, we caution 
that the enrichment scores are low, which precludes making 
confident conclusions (Table 5). 
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b) 



Suppressed 

Gefitinib Not Gef. 



Induced 

Gefitinib Not Gef. 



Etlotinib 




miR-7 



Figure 3. Venn diagrams of genes regulated by the kinase Inhibitors. 'Kinase Inhibitors' ovals represent genes regulated by all kinase 
inhibitors, 'Not Gef.' those by all inhibitors except Gefitinib. c) Venn diagrams of genes regulated by antibodies and by miR-7. 'All' represents the gens 
from all 346 samples, 'Kinase Inhibitors' ovals represent genes regulated by all kinase inhibitors. 
doi:1 0.1 371/journal.pone.01 02466.g003 



Effects of miR-7 

miR-7 suppresses EGFR production by direct interaction with 
two sites within the EGFR mRNA [24]. Pleiotropic, miR-7 
inhibits multiple targets downstream from and in addition to 
EGFR signaling, including Akt, ERKl/2, RAFl, IRS1&2 and 
PAKl. 

Predictably, the effects of miR-7 differ significantly from those 
of other agents. Ontological categories suppressed by miR-7 
notably include cell-cycle (Table 6). Ribosome assembly and 
biogenesis is, however, induced by miR-7 (Table 6). The 
MAPK3K family targets the miR-7 regulated genes; MAP3K14 
and RPK3 spc-cifically the miR-7 induced genes. However, 
because miR-7 inhibits multiple targets, it is currently not possible 
to assign its effects specifically to the EGFR-suppressing function. 

Discussion 

This work demonstrates the advantage of metaanalysis over 
single studies: metaanalysis provided 5 times more regulated genes 
than the largest single study. Importantly, coherent, single 
platform metaanalysis has advantages over an assortment of 
platforms, but in general analyses of large data sets provide more 
regulated genes than of smaller ones. Because we used free, 
publically available metaanalysis programs, this work can serve as 
a paradigm for integration and metaanalysis of transcriptional 
data in public repositories (Fig. 1). 

Large lists of regulated genes allowed us to identify novel 
ontological categories affected by EGFR inhibition. As expected, 
the suppressed genes are associated with cell-cycle, migration, 
transcription and protein synthesis, while the induced genes 
include ones associated with apoptosis, and inhibition of 



transcription and translation. Unexpectedly, the induced catego- 
ries also include genes associated with lysosome and with steroid 
hormone receptor activity. The induction of lysosomal genes by 
EGFR inhibition is a component of the autophagy, a process 
commonly associated with EGFR inhibition [25]. The induction 
of lysosomal genes may also play a significant role in, e.g., reducing 
the effectiveness of the inhibitors by degradation. 

Separate analyses of kinase inhibitors and antibodies identified 
important differences and commonalities. For example, antibodies 
suppress cell migration genes, much less the cell-cycle genes, while 
the reverse is true for kinase inhibitors. The differences do not 
derive from differences in cell types targeted, i.e., muscle or 
neuronal, they seem to be specific consequences of using different 
agents to inhibit EGFR. The molecular mechanisms causing these 
differences are not known; we speculate that antibodies, being 
large molecules, perturb the agglomeration of EGFR in the cell 
membrane [26]. This may affects the interaction between EGFR 
and other membrane-bound proteins transducing some of the 
downstream signals. 

Our results fit well with previous findings that combination 
treatments with multiple agents can have synergistic effects [13]. 
They suggest that combination treatments targeting EGFR, i.e., 
simultaneous use of antibodies and kinase inhibitors, may be 
beneficial for avoiding development of resistance. Moreover, we 
suggest that using specific combinations of agents can be fine- 
tuned and personalized to achieve patient-specific treatment 
responses. For example, we speculate that highly proliferative 
but rarely metastasizing cancers may benefit more from treatment 
with Erlotinib, which strongly affects cell cycle progression genes 
(Table 2), whereas highly metastatic tumors may benefit more 
from antibody therapies, which strongly affect cell motility 
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Table 6. Clusters of ontological categories regulated by miR-7. Note that the cell-cycle related clusters are suppressed by miR-7, 



while the transcriptional machinery is induced. 







Suppressed by miR-7 






By miR-7 only 






Enrichment Score: 8.97 


p-Val 




Enrichment Score: 2.48 


p-Val 


c 


mitotic cell cycle. 


3.E-1 1 




mitosis. 


3.E-03 




Enrichment Score: 8.06 




C 


M phase of mitotic cell cycle. 


3.E-03 


c 


M phase of mitotic cell cycle. 


1 .E-09 




Enrichment Score: 2.23 






Enrichment Score: 1.71 






centrosome. 


2.E-03 


N 


nuclear lumen. 


4.E-02 


c 


microtubule organizing center. 


4.E-03 




Enrichment Score: 1.56 






Enrichment Score: 1.83 




c 


mitotic sister chromatid segregation. 


6.E-03 




negative regulation of protein metabolism. 


7.E-03 




Enrichment Score: 1.42 




p 


regulation of protein metabolic process. 


8.E-02 




Heat shock protein DnaJ. 


2.E-02 




Enrichment Score: 1.37 












response to endogenous stimulus. 


1 .E-02 








s 


response to hormone stimulus. 


3.E-02 


Induced 




Enrichment Score: 16.22 


p-Val 




Enrichment Score: 17.26 


p-Val 


T 


structural constituent of ribosome. 


3.E-18 




translational elongation. 


3.E-24 




Enrichment Score: 15.82 




T 


Ribosome. 


1.E-18 


T 


cytosolic ribosome. 


4.E-18 




Enrichment Score: 15.21 






Enrichment Score: 3.35 






cytosolic ribosome. 


2.E-17 


T 


rRNA processing. 


1 .E-04 


T 


ribosomal subunit. 


9.E-16 




Enrichment Score: 2.61 






Enrichment Score: 3.75 






mitochondrial inner membrane. 


2.E-04 




rRNA processing. 


4.E-05 




Enrichment Score: 2.40 




T 


ribosome biogenesis. 


2.E-04 




Ubiquitin. 


5.E-03 




Enrichment Score: 2.99 






Enrichment Score: 1.94 






mitochondrial inner membrane. 


8.E-05 




intermediate filament cytoskeleton. 


7.E-02 




mitochondrial envelope. 


1 .E-03 




Enrichment Score: 1.93 






Enrichment Score: 2.59 






mitochondrial ATP synthase. 


5.E-04 




Ubiquitin conserved site. 


1 .E-03 




Enrichment Score: 1.78 




P 


Ubiquitin. 


3.E-03 




mitochondrial lumen. 


l.E-01 




Enrichment Score: 2.13 












mitochondrial ATP synthase complex. 


3.E-04 










proton-transporting ATP synthase. 


5.E-04 
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(Table 5). Gefitinib seems less proapoptotic than other kinase 
inhibitors, which may reduce side effects of targeting EGFR in 
specific cases (Table 3). These are just guidelines, which will need 
experimental corroboration. 

Antibodies specifically suppress the developmental effects of 
EGFR (e.g., in skin, hair, vasculature and gonads). Antibodies also 
specifically induce genes associated with the contractile apparatus. 
Such effects perhaps depend on altered EGFR-containing multi- 
protein complex formation in cell membrane, a cytoskeleton- 
dependent process [4,26]. 

While kinase inhibitor drugs are supposed to act with same 
mechanisms, there are characteristic differences in on-off rates, 
receptor conformation and accessibilities of its serine/ threonine/ 
tyrosine target substrates [27]. Apparently, individual inhibitors 
use distinct mechanisms to achieve similar results: cell-cycle 
inhibitors are induced by Gefitinib, conversely cell-cycle promot- 
ers are suppressed by Erlotinib - the same results achieved by 



different routes. Targets of GSK3 are significantly suppressed by 
Gefitinib, not by other agents. Neural tube closure and steroid 
hormone receptor activity are particular targets of induction by 
Gefitinib, which deserves further studies. These results may 
suggest specific preference for use of Gefitinib in certain tumors, 
e.g., glioblastomas, and more problematic side-effects in other 
tumors. 

The signal transducing kinases inhibited by EGFR-targeting 
agents largely overlap. They include known members of the 
EGFR signaling cascade, MAP3K being the most prominent. 
Although individual agents are associated with individual kinases, 
the differences seem subtle and it is unclear presentiy whether this 
is due to statistical effects of different study sizes, i.e., numbers of 
samples, or real mechanistic differences responding to different 
agents. Even more homogenous are the transcription factors that, 
presumably, control the expression of the regulated genes: similar 
sets of TFs respond to multiple EGFR inhibitors, although there 
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are difTerences in their relative orders and p-values. The 
significance of these subde differences is at present dubious. 
Highly overlapping sets of TFs regulate both the suppressed and 
the induced genes. This would suggest that constellations of TFs in 
the promoters of regulated genes, their interactions with auxiliary 
proteins or as yet unknown TFs determine whether a gene is 
induced or suppressed by the inhibition of EGFR. 

Supporting Information 

Table SI The 2537 genes suppressed by the EGFR 
inhibitors. The Afiymetrix ID, gene name, cytoband, Entrez ID, 
associated KEGG pathway, official symbol and OMIM disease 
associations are included. 
(XLSX) 

Table S2 The 2251 genes induced by the EGFR inhib- 
itors. 

PCLSX) 

Table S3 Ontological categories. Top 100 ontological 
categories, by p-values, are listed for each list of genes mentioned 
in the text. 
pCLSX) 

Table S4 Kinases associated with regulated genes. The 

MAPK3 family kinases participate in both induction and 
suppression of genes by EGFR inhibitors; they are particularly 
important in the induction by miR-7. Note also the differences 
between the kinases responsive to Gefitinib and those responsive to 

References 

1. Schneider MR, Werner S, Paus R, Wolf E (2008) Beyond wavy hairs: the 
epidermal growth factor receptor and its ligands in skin biology and pathology. 
AmJ Pathol 173: 14-24. doi: 10.2353/ajpath.200a.070942. 

2. Bae JH, Schlessingrr J (2010) A.symmetric tyrosine kinase arrangements in 
activation or autophosphorvlation of receptor tvrosine kinases. Mol Cells 29: 
443-448. doi: 410.1007/sl00,i9-10010-10080-100.5,x 

3. Ullrich A, Schlessinger J (1990) Signal transduction by receptors with tvrosine 
kinase activity. Cell 61: 203-212. 

4. Chung I, Akita R, Vandlen R, Toomre D, Schlessinger J, ct al. (2010) Spatial 
control of EGF receptor activation by reversible dimerization on living cells. 
Nature 464; 783-787. doi: 710.1038/nature08827. 

5. Wheeler DL, Dunn EF, Harari PM (2010) Understanding resistance to EGFR 
inhibitors-impact on future treatment strategies. Nat Rev Clin Oncol 7: 493- 
507. 

6. Singh AB, hlarris RC (200,^) Autocrine, paracrine and juxtacrine signaling by 
EGFR ligands. Cell Signal 17: 1183-1193. 

7. Takeuchi K, Ito F (2010) EGF receptor in relation to tumor development: 
molecular basis of responsiveness of cancer cells to EGFR-targeting tyrosine 
kinase inhibitors. Febs J 277: 316-326. 

8. Zheng Y, Zhang C, Croucher DR, Soliman MA, St-Denis N, et al. (2013) 
Temporal regulation of EGF signalling networks bv the scaffold protein Shcl. 
.Nature 499: 166-171. doi: 1 10. 1038/naturel2308. ' 

9. Rull-Jamison S, Chen K, Cohen S (1993) Induction by EGF and interferon- 
gamma of tyrosine phosphorylated DNA binding proteins in mouse liver nuclei. 
Science 261: 1733-1736. 

10. Pao W, Chmielecki J (2010) Rational, biologically based treatment of EGFR- 
mutant non-small-cell lung cancer. Nat Rev Cancer 10: 760-774. doi: 
710.103B/nrc2947. 

1 1 . Pomerantz RG, Mirvish ED, Geskin LJ (201 0) Cutaneous reactions to epidermal 
growth factor receptor inhibitors. J Drugs Dermatol 9: 1229—1234. 

12. Ciardiello F, I'ortora G (2008) EGFR antagonists in cancer treatment. 
N EnglJ .\Ied 35R: 1160-1174. doi: 1 1 10. 10,i6/NEJ.\Ira0707704. 

13. Perera RM, Narita Y, Furnari FB, Gan HK, Murone C, et al. (200.')) I'reatment 
of human tumor xenografts with monoclonal antibody 806 in combiiiation with 
a prototypical epidermEd growth factor receptor-specific antibody generates 
enhanced antitumor activity. Clin Cancer Res 11: 6390-6399. 

14. Hong F, BreitUng R, McEntee CW, Wittaer BS, Nemhauser JL, ct al. (2006) 
RankProd: a bioconductor package for detecting diiferentially expressed genes in 
meta-analysis. Bioinformaties 22: 282.')-2827. 

15. Mimoso C, Lee DD, Za\'adilJ, Fomic-Canic Al, Blumenberg M (2013) Analysis 
and Meta-analysis of Transcriptional Profiling in Human Epidermis. Methods 
Mol Biol 3: 3. 



Erlotinib, such as GSK3s. No kinase reached statistical signifi- 
cance in the analysis of genes regulated by antibodies. 

(XLSX) 

Table S5 Clusters of ontology categories specifically 
regulated by individual kinase inhibitors. Note the much 
relaxed enrichment scores used as cut-offs because of relatively 
short lists of genes. 
(XLSX) 

Table S6 Top 30 transcription factors associated with 
the regulated genes. Grey marks the top 10 factors in analysis 
of all 346 samples; the same 10 are marked with grey in the 
individual analyses as well, a) Suppressed genes, b) Induced genes. 
We find that 8 out of 10 top transcription factors regulate both the 
suppressed and the induced genes, these are marked with asterisks. 
The transcription factors associated with the miR-7 regulated 
genes and with genes suppressed by the antibodies did not reach 
statistical significance. 
(XLSX) 

Text SI Description of studies used in this metaana- 
lysis. 

(DOCX) 

Author Contributions 

Conceived and designed the experiments: MB. Performed the experiments: 
MB. Analyzed the data: MB. Contributed reagents/materials/analysis 
tools: MB. Contributed to the writing of the manuscript: MB. 



16. Gautier L, Cope L, Bolstad BM, Irizarry RA (2004) affy-analysis of Afiymetrix 
GeneChip data at the probe level. Bioinformaties 20: 307-315. 

17. Lachmann A, Ma'ayan A (2010) Lists2Networks: integrated analysis of gene/ 
protein lists. BMC Bioinformaties 11: 87. 

18. Dennis G, Jr, Sherman BT, Hosack DA, Yang J, Gao W, et al. (2003) DAVID: 
Database for Annotation, Visualization, and Integrated Discovery. Genome Biol 
4: P3. 

19. Liu V, M.rjumdrr S, ,\leCall \V, Sartor CI. .Mohlrr JL, el al. f200:)) Inhibnion of 
HER-2/neii kinase impairs androgen receptor recruitment to the androgen 
responsi\'e enhancer. Cancer Res 65: 3404-3409. 

20. Stabile LP, Lyker JS, Gubish CT, Zhang W, Grandis JR, et al. (2005) Combined 
targeting of the estrogen receptor and the epidermal growth factor receptor in 
non-small cell lung cancer shows enhanced antiproliferative effects. Cancer Res 
65: 1459-1470. 

21. Kuo LC, Lin PC, Wang KF, Yuan MK, Chang SC (201 1) Successful ti-eatinent 
of gefitinib-induced acute interstitial pneumonitis with high-dose corticosteroid: 
a case report and literature review. Med Oncol 28: 79-82. doi: 10.1007/sl2032- 
12010-19424-12034. 

22. Ali A, Hoeflich KP, Woodgett JR (200 1) Glycogen synthase kinase-3: properties, 
functions, and regulation, (^hem Re\- 101: 2527-2540. 

23. Casas-Selves .\1. KinrJ. Zhani; Z, Helfrich BA, Gao D, et al. (2012) Tankyrase 
and the canonical Wnt pathway protect lung cancer cells from EGFR inhibition. 
Cancer Res 72: 4154-4164. doi: 41 10. 1 158/0008-5472. CAN-41 1 1-2848. 

24. Kalinowski FC, Giles KM, Candy PA, Ali A, Ganda C, et al. (2012) Regulation 
of epidermal growth factor receptor signaling and erlotinib sensitivity in head 
and neck cancer ccUs by miR-7. PLoS One 7: e47067. doi: 47010.41371/ 
journal.pone.0047067. 

25. Dragowska WH, Weppler SA, WangJC, Wong LY, Kapanen Al, et al. (2013) 
Induction of autophagy is an early response to gefitinib and a potential 
therapeutic target in breast cancer. PLoS One 8: e76503. doi: 76510.71371/ 
journal.pone.0076503. eCollection 0072013. 

26. Clayton AH, Tavarnesi ML, Johns IXi (2007) Unligated epidermal growth 
factor receptor forms higher order oligomers within microclusters on A43 1 cells 
that are sensitive to tyrosine kinase inhibitor binding. Biochemistry 46: 4589— 
4597. 

27. Wood ER, Truesdale AT, McDonald OB, Yuan D, HasseU A, et al. (2004) A 
imique structure for epidermal grovrth factor receptor bound to GW572016 
(Lapatinib): relationships among protein conformation, inhibitor off-rate, and 
receptor activity in tumor cells. Cancer Res 64: 6652-6659. 



PLOS ONE I www.plosone.org 



14 



September 2014 | Volume 9 | Issue 9 | e102466 



